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ABSTRACT 

In order to understand the heating of the solar corona it is crucial to obtain empirical information on the 
magnetic field in its lower boundary (the transition region). To this end, we need to measure and model the 
linear polarization produced by scattering processes in strong UV lines, such as the hydrogen Lya line. The 
interpretation of the observed Stokes profiles will require taking into account that the outer solar atmosphere 
is highly structured and dynamic, and that the height of the transition region may well vary from one place in 
the atmosphere to another. Here we report on the Lya scattering polarization signals we have calculated in a 
realistic model of an enhanced network region, resulting from a state-of-the-art radiation MHD simulation. This 
model is characterized by spatially complex variations of the physical quantities at transition region heights. 
The results of our investigation lead us to emphasize that scattering processes in the upper solar chromosphere 
should indeed produce measurable linear polarization in Lya. More importantly, we show that via the Hanle 
effect the model's magnetic field produces significant changes in the emergent Q/I and U/I profiles. Therefore, 
we argue that by measuring the polarization signals produced by scattering processes and the Hanle effect in 
Lya and contrasting them with those computed in increasingly realistic atmospheric models, we should be able 
to decipher the magnetic, thermal and dynamic structure of the upper chromosphere and transition region of 
the Sun. 

Subject headings: polarization — scattering — radiative transfer — Sun: chromosphere — Sun: transition 
region — Sun: surface magnetism 



1. INTRODUCTION 

In order to understand the heating of the solar corona it is 
crucial to obtain empirical information on the strength and 
orientation of the magnetic field in its lower boundary (the 
transition region), where over very short distances (~100 km) 
the temperature suddenly rises from < 10 4 K to > 10 6 K and 
the atmospheric pl asma changes from partia lly to practically 
fully ionized (e.g jGolub & Pasach off 2010). To this end, we 
need to measure and model observables sensitive to the mag- 
netic field of the sola r transition region. 

In a recent paper, Trujillo Buen o et al.l (1201 ll) argued that 
the hydrogen Lya line is expected to show measurable scat- 
tering polarization when observing the solar disk, and that 
via the Hanle effect the line-center linear polarization ampli- 
tude shows a good sensitivity to magnetic field strengths be- 
tween 10 and 100 G, approximately. These conclusions were 
reached through detailed radiative transfer calculations in 
semi-empirical and hydrodynamical models of the solar atmo- 
sphere, assuming complete frequency redistribution (CRD) 
and neglecting quantum interference between the two upper 
levels of the Lya line (see Step an & Trujillo Buenold201 lb for 
details on the atomi c model and numeri cal method of solu- 
tion). As shown by Bell uzzi et al.l (1201 2l) . these last two ap- 
proximations are suitable for estimating the polarization sig- 
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nals at the core of Lya, which is the spectral line region where 
the Hanle effect in Lya operates. 8 

The above-mentioned investigations were carried out using 
one-dimensional models of the exten ded solar atmo s phere , 
such as the semi-empirical models of iFontenla et alJ (1 19931) 
and the hydrodynamical models of Carlsso n & Steinl d 1997b . 
However, both o bservations (e.g., IVourlidas et al.l [20 1 Oh and 
simulations ("e.g.. iLeenaarts et al1l2012l) show that the outer 
solar atmosphere is highly structured and dynamic, departing 
radically from a uniform, plane-parallel configuration. Thus, 
the height of the transition region must vary from one place 
in the atmosphere to another, most probably delineating a 
highly corrugated surface. It is in this type of spatially com- 
plex plasmas where strong transition region lines like Lya 
originate. The presence of horizontal atmospheric inhomo- 
geneities in a stellar atmosphere breaks the axial symmetry of 
the radiation field within the medium, and this type of sym- 
metry breaking may produce changes in the scattering polar- 
ization signals that co uld be confused with those caused by 
the Hanle effect (e.g., [M anso Sainz & Trujillo Bueno 1201 lb 



lAnusha & N agendral l201 ll) . It is therefore important to in 
vestigate the problem of scattering polarization and the Hanle 
effect in Lya using spatially complex models of the outer solar 
atmosphere, where the physical quantities may vary abruptly 
not only along the radial direction but also horizontally. To 
this end, Stepan & T ruji llo Buenol d2012 a b) have developed 
a three-dimensional (3D), multilevel radiative transfer code 
for doing simulations of the spectral line polarization caused 
by scattering processes and the Hanle and Zeeman effects 
within the f ramework of the quantum theory of pola rization 
described in Landi DeglTnnocenti & Landolfi (2004). 
The present paper represents a first step towards that goal. 

8 Partial frequency redistribution (PRD) and /-state interference can, how- 
ever, produce large linear polarization signals in the wings o f the Q/I pro- 



file o f Lyo-, especially for close-to-the-limb line of sights (see Belluzzi et al. 
|20H . 
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Fig. 1 . — The horizontal variation of the kinetic temperature of the 3D model 
at a height of 2000 km. The dashed line indicates the X-position of the Z-Y 
slice chosen for obtaining the 2D model studied in this Letter. 

Given the complexity of the radiative transfer problem of the 
Hanle effect in 3D model atmospheres here we start by con- 
sidering a two-dimensional (2D) model atmosphere character- 
ized by spatially complex variations of the physical quantities 
at chromospheric and transition region heights, taken from 
the 3D radiation magneto-hydrodynamic (MHD) simulations 
mentioned below. Of particular interest here is to compare 
the Q/I and U/I profiles computed assuming B = G at each 
spatial grid point of the considered MHD atmospheric model 
with the fractional linear polarization signals calculated tak- 
ing into account the Hanle effect caused by the model's mag- 
netic field. 

2. FORMULATION OF THE RADIATIVE TRANSFER PROBLEM 

We have selected a snapshot from a state-of-the-art ra- 
diation MHD simula tion performed with the Bifrost code 
( Gudiksen et alj 1201 lb taking into account non-equilibrium 
hydrogen ionization. A snapshot from the same simulation 
(but 30 seconds earlier in time), which has a magnetic field 
configuratio n representative of an enhanced network region, 
was used bv lLeenaarts et alj d2012l) and we refer the reader to 
that paper for details on the 3D structure. Figure Q] illustrates 
the snapshot we have selected from this simulation. The mag- 
netic field has a predominantly bipolar structure, with mag- 
netic field lines connecting two clusters of photospheric mag- 
netic concentrations of opposite polarity and reaching chro- 
mospheric and coronal heights. In this Letter we consider a 
two-dimensional (2D) model atmosphere extracted from this 
3D snapshot; its physical quantities are those found within the 
Z-Y plane at X=12 Mm (see the dashed line in Fig.Q]). Since 
the magnetic field lines in this 3D snapshot are approximately 
perpendicular to such a Z-Y plane and the temperature be- 
tween the magnetic patches shows filamentary structure that 
follows the general orientation of the magnetic field, we be- 
lieve that the 2D approximation is reasonable for computing 
the emergent Stokes profiles for line of sights (LOS) with az- 
imuths^- = 0° or^- = 180° (i.e., for LOS contained in any 
Z-X plane). 

The Lya line results from two blended transitions between 
the ground level of hydrogen, ls 2 Si/2, and the 2p 2 Pi/2 and 
2p 2 P3/2 excited levels. The only level that contributes to 



the scattering polarization and the Hanle effect of Lya is the 
2p 2 f > 3/2 (e.g.. iTrujillo Bueno et alJl201 ll). We hav e adop ted 
the same atomic data used by iTruiillo Bueno et alj d201 ll) in 
their one-dimensional (ID) investigation, the only difference 
being that we have used a model atom composed of the above- 
mentioned levels and the 2s 2 Si/2 level whose collisional cou- 
pling with the P-levels may contribute to a slight depolar- 
ization of the line. As shown by Ste pan & Truiillo Buenol 
( 12011b this is a very suitable atomic model for computing 
the line-center polarization of Lya. We aim at comparing 
the Q/I and U/I profiles computed assuming B = G at 
each spatial grid point of the above-mentioned 2D atmo- 
spheric model with the fractional linear polarization signals 
calculated taking into account the Hanle effect caused by 
the model's magnetic field. To this end, we have applied 
the radiati ve transfer code mentioned in Se ction 1. As ad- 
vanced by Step an & Trujillo Buenol (1201 2ab . two novel fea- 
tures of this 3D code for the transfer of polarized radia- 
tion is a formal solver based on short-characteristics with 
quadratic Bezier-splines interpolation (which implies a signif- 
icantly more stable formal solver than those based on the stan- 
dard parabolic interpolation method) and a highly-convergent 
iterative m ethod based on the non- l inear multigrid method de- 
scribed in lFabiani Bendicho et alJ (1997). For the angular in- 
tegration we have used Gaussian quadrature for the ray incli- 
nations with respect to the vertical Z-axis (with 7 inclinations 
per octant) and the trapezoidal rule for the ray azimuth (with 
5 azimuths per octant). 

The top panel of Fig. [2] shows the spatial variation of 
the kinetic temperature in the chosen 2D model atmosphere, 
where it can be clearly seen that the height of the transi- 
tion region changes drastically from one horizontal location 
to another, outlining a highly crumpled surface. The next 
panel of Fig. [2] shows the spatial variation of the fractional 
anisot ropy of the Lya radiation, M /J% (see equations (2) and 
(3) in iTrujillo Bueno et al.ll2.01 lb . calculated after obtaining 
with the above-mentioned radiative transfer code the self- 
consistent solution to the problem of resonance line polar- 
ization and the Hanle effect. As shown in this panel, the 
height variation of J%/J%, at most horizontal positions Y, is 
qualitatively similar to tha t shown for Lyes' in figure 1 of 
Truiill oBueno et all (1201 lb . which corresponds to the solu- 
tion of the same non-LTE polarization transfer problem in the 
semi-empirical model C of Fontenla et al. (1993; hereafter, 
FAL-C model). Note that at most horizontal positions the 
fractional anisotropy is zero at all photospheric and chromo- 
spheric heights, but that it suddenly becomes significant right 
at the location of the transition region, with larger negative 
values than in the case of the semi-empirical FAL-C model 
atmosphere. This is important because of the following two 
reasons: (1) the corrugated surface that delineates the model's 
transition region practically coincides with that correspond- 
ing to line-center optical depth unity along the line of sight 
(LOS), and (2) the line-center amplitude of the Q/I profile 
of the Lya line is approximately proportional to the J^/J^ 
value at such transition r egion heights (cf., equation (1) of 
ITruiillo Bueno et alj|2011b . 

The extra three panels of Fig.|2]show the spatial variation of 
the strength, inclination and azimuth of the model's magnetic 
field. Note that at the crinkled surface that outlines the lo- 
cation of the model's transition region (which approximately 
coincides with the height where the Lya line-center optical 
depth is unity along the LOS) the strength of the magnetic 
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Fig. 2. — Physical quantities across the vertical 2D slice denned in Fig.fT] From top to bottom: (a) temperature, (b) the fractional anisotropy of the hya radiation 
field (the black contours show the plasma ft parameter), (c) magnetic field strength B, (d) inclination of magnetic field with respect to the positive Z-axis, and 
(e) azimuth of magnetic field with respect to the positive X-axis. The white line in each panel indicates the height, associated to each horizontal Y-point, where 
the line-center optical depth is unity along a LOS with /j = 0.3, and we point out that it is very similar to that corresponding to other /i values. Note that in 
one-dimensional semi-empirical models the transition region to the coronal temperatures occurs instead at a rather precise height (e.g., at about 2200 km in the 
FAL-C model), which approximately coincides with the height where the Lya line-center optical depth is unity along the LOS. In both types of atmospheric 
models (MHD and semi-empirical) the line-core hya radiation originates around such transition region heights. 



field varies between a few gauss and 60 gauss, approximately. 
Since the Hanle effect in Lya is mainly sensitive to magnetic 
streng ths between 10 and 100 G (cf., iTrujillo Bueno et all 
1201 ll) we should expect a significant impact on the emergent 
Q/I and U/I line-core signals. 

3. THE EMERGENT Q/I AND U/I SIGNALS 

The first step needed to compute the emergent Q/I and U/I 
signals is to obtain, at each spatial point within the chosen 
model atmosphere, the self-consistent values of the density- 



matrix elements corresponding to each j-level of our 4-level 
atomic model (i.e., 9 pg elements quantifying the overall pop- 
ulation of each j'-level and the population imbalances and co- 
herences among the sublevels of the 2p 2 Pi/2 level). To this 
end, we have applied the 3D radiative transfer code developed 
by Ste pan & Tru iillo Bueno] d2012alfbl) for solving the two fol- 
lowing cases: (a) the zero-field reference case, by forcing 
B — G at each spatial grid point in the atmospheric model 
(top panels) and (b) the Hanle effect case, by taking into ac- 
count the model's magnetic field. The general multilevel ex- 
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Fig. 3. — The emergent Q/I and U/I profiles calculated at each Y-point in 
the 2D model of Fig. [2] for a LOS with fi = 0.3 and x = 0°. Top panels: 
assuming B = G at each point within the 2D model. Bottom panels: taking 
into account the Hanle effect produced by the model's magnetic field. The 
Stokes profiles at any given Y position were averaged over a 5" interval. 

pressions of the radiative transfer coefficients and of the sta- 
tistical equilibrium equations corresponding to this so-called 
N LTE problem of the second kind can be fo und in Section 7.2 
of lLandi DeglTnnocenfi & Landolfil d2004l) . 

In Fig. [3] we show, for a LOS with fi - 0.3 and azimuth 
X = 0°, the emergent Q/I and U/I profiles calculated for each 
of the above-mentioned two cases. The first important point 
to emphasize is that the amplitudes of the linear polarization 
signals are very significant (~1%), both in the absence (top 
panels) and in the presence (bottom panels) of the model's 
magnetic field, and that they show significant spatial varia- 
tions. 

Note that Q/I is always negative and that its amplitude 
is significantly larger in the zero field case. The fact that 
U/I + for the B — 0G case (top panels) indicates that the 
model's horizontal inhomogeneities (e.g., in temperature and 
density) are efficient in breaking the axial symmetry of the 
incident radiation field at each point within the medium. In- 
terestingly enough, both the amplitude and the sign of U/I are 
significantly modified by the action of the Hanle effect (com- 
pare the top and bottom panels), mainly in the central part 
of the 2D model where the magnetic field strength is signifi- 
cantly larger than in the model's boundaries (see Fig. |2j. 

4. THE MAGNETIC SENSITIVITY OF THE SPATIALLY-AVERAGED 
Q/7 AND U/I SIGNALS 

It is also interesting to consider the spatially-averaged 
Q/I and U/I signals, including their sensitivity to the mean 
magnetic field strength obtained by averaging the magnetic 
strength at each point of the corrugated surface that defines the 
spatial location of the model's transition region. To this end, 
we have multiplied the strength of the local magnetic field by 
a scaling factor < / < 8. Figure H] shows the results for an 
LOS with n = 0.3. 

The dashed red line of Fig. |4] shows the magnetic sen- 
sitivity of the total fractional linear polarization, Pl = 

V<fi> 2 + (U) 2 /(I), where it can be seen that the model's mag- 
netic field produces via the Hanle effect a significant depo- 
larization with respect to the zero-field reference case. More 
interesting is to note that (U)/(I) ~ in the absence of mag- 
netic fields while (U)/{I) ~ +1% with the model's magnetic 
field (see the black dotted lines). As shown by the blue solid 
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Fig. 4.— The <£>/</>, (U)/(I) and P L = vW + (U) 2 /(I) spatially- 
averaged amplitudes versus the mean field strength of the model's transition 
region, as discussed in Section 4. The LOS chosen have /.( = 0.3 and az- 
imuths^ = 0° (black lines) and^- = 180° (blue lines). The vertical green 
line indicates the mean field strength of the model's transition region. Other 
solutions have been obtained by scaling the modulus of the model's mag- 
netic field, with a scaling factor < / < 8. Note that the (Q)/(I) ampli- 
tude corresponding to the zero-field reference case in the chosen 2D model 
is significantly larger than that found by Trujillo Bueno et al. (2011) in the 
one-dimensional FAL-C model. 

lines, the Hanle effect of the model's magnetic field has a sig- 
nificant impact on the (Q)/(I) amplitude, whose sign remains 
the same for the two considered LOS azimuths. 

5. CONCLUDING COMMENTS 

In order to obtain quantitative information on the magnetic 
field of the upper chromosphere and transition region of the 
Sun we need to measure and model the linear polarization 
caused by scattering processes in some strong resonance lines, 
such as the hydrogen Lye line at 1216 A. This is because via 
the Hanle effect such polarization is sensitive to the strength 
and orientation of the magnetic field expected for the solar 
transition region plasma (with B < 100 G). The interpreta- 
tion of the observed Stokes profiles will require taking into 
account that the outer solar atmosphere is highly structured 
and dynamic. For this reason it is important to do numerical 
simulations of the scattering polarization that is produced in 
realistic MHD models of the solar atmosphere, paying par- 
ticular attention to understanding the impact of the model's 
magnetic field on the emergent Q/I and {///profiles. 

In this Letter we have made a first step towards that goal, by 
solving the multilevel radiative transfer problem of the scat- 
tering polarization and Hanle effect in Lya in a 2D snapshot 
taken from the 3D radiation MHD simulation of an enhanced 
network region described in Leenaarts et al. (2012). To this 
end, we applied the 3D radiative transfer code described in 
IStepan & Trujillo Buenol (l2012albh . 

Our model calculations of the line-core Q/I and U/I sig- 
nals at each point on the model's surface show that the hori- 
zontal inhomogeneities of the chosen atmospheric model pro- 
duce a significant impact on the local scattering polarization 
amplitudes, and that these linear polarization signals are sig- 
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nificantly different in the absence and in the presence of the 
Hanle effect produced by the model's magnetic field. More- 
over, the spatially-averaged Q/I and U/I signals show also 
an interesting sensitivity to the model's magnetic field. We 
therefore conclude that by measuring the polarization signals 
produced by scattering processes and the Hanle effect in Lya 
and contrasting them with those computed in increasingly re- 
alistic atmospheric models, we should be able to decipher the 
magnetic, thermal and dynamic structure of the upper chro- 
mosphere and transition region of the Sun. 

We are grateful to Rafael Manso Sainz (IAC) for carefully 
reviewing this Letter and for several interesting discussions. 
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